KILLED BROWNIAN MOTION WITH A PRESCRIBED 
LIFETIME DISTRIBUTION AND MODELS OF DEFAULT 



BORIS ETTINGER, STEVEN N. EVANS, AND ALEXANDRU HENING 

Abstract. The inverse first passage time problem asks whether, for a Brow- 
nian motion B and a nonnegative random variable there exists a time- 
varying barrier b such that ¥{Bs > b(s), < s < t} = P{C > t}. We study 
a "smoothed" version of this problem and ask whether there is a "barrier" b 
such that E[exp(-A /g ip{Bs - b{s)) ds)] = ¥{( > t}, where A is a killing rate 
parameter and i/) : R — > [0, 1] is a non-increasing function. We prove that if tp 
is suitably smooth, the function t i— > ¥{( > t} is twice continuously differen- 
tiable, and the condition < _£Ll££Mi>ii < \ holds for the hazard rate of 
then there exists a unique continuously differentiable function b solving the 
smoothed problem. We show how this result leads to flexible models of default 
for which it is possible to compute expected values of contingent claims. 



1. Introduction 

Investors are exposed to credit risk, or counterparty risk, due to the possibility 
that one or more counterparties in a financial agreement will default; that is, not 
honor their obligations to make certain payments. Counterparty risk has to be 
taken into account when pricing a transaction or portfolio, and it is necessary to 
model the occurrence of default jointly with the behavior of asset values. 

The default time is sometimes modeled as the first passage time of a credit index 
process below a barrier. Black and Cox jBC76| were among the first to use this 
approach. They define the time of default as the first time the ratio of the value of 
a firm and the value of its debt falls below a constant level, and they model debt 
as a zero-coupon bond and the value of the firm as a geometric Brownian motion. 
In this case, the default time has the distribution of the first-passage time of a 
Brownian motion (with constant drift) below a certain barrier. 

Hull and White [HWOlj model the default time as the first time a Brownian 
motion hits a given time-dependent barrier. They show that this model gives the 
correct market credit default swap and bond prices if the time-dependent barrier 
is chosen so that the first passage time of the Brownian motion has a certain 
distribution derived from those prices. Given a distribution for the default time, it 
is usually impossible to find a closed-form expression for the corresponding time- 
dependent barrier, and numerical methods have to be used. 



Date: June 27, 2012. 

2010 Mathematics Subject Classification. 60J70, 91G40, 91G80. 

Key words and phrases, credit risk, inverse first passage time problem, killed Brownian motion. 
Cox process, stochastic intensity, Feynman-Kac formula. 
SNE supported in part by NSF grant DMS-0907630. 



2 



BORIS ETTINGER, STEVEN N. EVANS, AND ALEXANDRU HENING 



We adopt a perspective similar to that of |HW01j . Namely, we model the default 
time as 

/■* 

(1.1) 



T := inf a > : A / il){Ys - b{s)) ds > U 



where the diffusion Y is some credit index process, U is an independent mean one 
exponentially distributed random variable, < 5; 1 is a suitably smooth, non- 
increasing function with limj;_>._oo i'ix) ~ 1 and limj;^+oo ^{x) ~ 0, and A > is 
a rate parameter. Then, 



(1.2) 



»{t ><} = E 



exp -A / ^{Ys-b{s))ds 



The random time r is a "smoothed-out" version of the stopping time of Hull and 
White; instead of of killing Y as soon at it crosses some sharp, time-dependent 
boundary, we kill Y at rate Xip{y — b{t)) if it is in state y G R at time t > 0. That 
is, 

Urn P{t G (t,t + At) I {Ys)o<s<t, T > t}/At = \ij{Yt - b{t)). 

When the credit index value Yt is large, corresponding to a time t when the coun- 
terparty is in sound financial health, the killing rate X4>{Yt — 6(t)) is close to and 
default in an ensuing short period of time is unlikely, whereas the killing rate is 
close to its maximum possible value. A, when Yt is low and default is more probable. 
Note that if we consider a family of [0, l]-valued, non-increasing functions tp that 
converges to the indicator function of the set {a; G K : a; < 0} and A tends to oo, 
then the corresponding stopping time r converges to the Hull and White stopping 
time inf{t > : F* < b{t)}. 

The hazard rate of the random time t is 



P{t e I r > t} 



dt 



(1.3) 



lim ■ 

At^O 



lim ■ 

At^O 



lim ■ 

AtiO 



'{r e {t,t + At)} 
AiP{T > t} 

' {a/o ib{Ys - b{s))ds <U< A/o+^* V(n - b{s)) ds} 
AtV {a /q ij{Ys - b{s)) ds < L/} 



Atl 



exp 



(-A/o*^(y.-5(s)), 



AE 


^(Ft -5(t))exp 


(-A/J^(n-6(s))ds)" 


E 


exp 


(-A/>(n-5(s))ds) 





On the other hand, suppose that C is a non-negative random variable with sur- 
vival function t H> G{t) := > t}. Writing g for the derivative of G, the 
corresponding hazard rate is 

As a result, a necessary condition for a function b to exist such that the correspond- 
ing random time r has the same distribution as ^ is that 

(1.4) < ~g{t) < XG{t),t>0. 
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We show in Theorem 12.11 that if F is a Brownian motion with a given suitable 
random initial condition, the assumption ()1.4p holds, and the survival function G 
is twice continuously differentiable, then there is a unique differentiable function 
h such that the stopping time r has the same distribution as C- In particular, we 
establish that the function h can be determined by solving a system consisting of 
a parabolic linear PDE with coefficients depending on h and a non-linear ODE for 
h with coefficients depending on the solution of the PDE. Note from (jl.2|) that 
changing the function on a set with Lebesgue measure zero does not affect the 
distribution of t, and so we have to be careful when we talk about the uniqueness 
of h. This minor annoyance does not appear if we restrict to continuous b. 

In Theorem 14.11 we give an analogue of the existence part of the above result 
when -0 is the indicator of the set {a;eM:a;<0}. 

Having proven the existence and uniqueness of a barrier h, we consider the pricing 
of certain contingent claims in Section [5] For simplicity, we take the asset price 
{Xt)t>Q to be a geometric Brownian motion 

—— = fidt + adWt , 
Xt 

where is a standard Brownian motion. We take the credit index {Yt)t>o to be 
given by 

dYt = dBt 

where B is another standard Brownian motion, and take the default time to be 
given by (jl.ip . where the exponential random variable U is independent of the 
asset price X and the credit index Y. We assume that the Brownian motions 
W and B are correlated; that is, that their covariation is [i3,W^]t = pt for some 
constant p g [—1, 1]. We consider claims with a payoff of the form F{Xt)1{t > T} 
for some fixed maturity T. We show how it is possible to compute conditional 
expected values such as 

F{Xt)1{t > T} I iX,)o<s<t, T > t 

In Section [6] we report the results of some experiments where we solved the 
PDE/ODE system for the barrier b numerically. Lastly, in Section [71 we follow 
[DPllj to demonstrate how it is possible to use market data on credit default swap 
prices to determine the survival function G. 

1.1. The FPT and IFPT problems. We end the this introduction with a brief 
discussion of the literature dealing with first passage times of diffusions across time- 
dependent barriers. 

Consider a Brownian motion {Bt)t>o defined on a filtered probability space 
{il,P,J-, (J^t)t>o) which satisfies the usual conditions. Define the diffusion {Yt)t>o 
via the SDE 

dYt^^I{Yt,t) dt + a{Yt,t) dBt, 

where we assume that the coefficients /i : M x IR+ M. and ct : M x IR+ M.+ are 
such that the SDE has a unique strong solution. 

For a Borel function b : M+ — )■ R := R U {±oo}, the first passage time of the 
diffusion process Y below the barrier b is the stopping time 

(1.5) f = inf{i >0:Yt< 6(t)}. 
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The following two problems related to this notion have been discussed in the liter- 
ature. 

The First Passage Time problem (FPT): For a given barrier h : R_|_ — > R, com- 
pute the survival function G of the first time that X goes below b; that is, find 

(1.6) G(t) P{f > i}, t>0. 

The Inverse First Passage Time problem (IFPT): For a given survival function 
G, does there exist a barrier b such that G{t) = P{f > t} for all t > 0? 

A large class of first passage time problems may be solved within a PDE frame- 
work. Let u{x,t) = -^FlYt < x,T > t} he the sub-probability density of the 
diffusion Y killed at f. Then, by the Kolmogorov forward equation, u satisfies 

{ut{x,t) = i(cr^u)j;^ - i^iu)^, X > b{t), t>0, 
u{x,t)=0, x<b{t),t>0, 
u{x,0)^f{x), xeR, 

where / is the probability density of Yg. For nice enough functions b this system 
has a unique solution and we can find the survival probability 

/•oo 

G{t)^F{f>t}^ u{x,t)dx, t>0. 
Jb{t) 

This approach is used in |Ler861 IValOQj to get closed form solutions for some classes 
of boundaries. An integral equation technique is used in jPes021 IPS06[ IValOQ] to 
find the derivative g{t) = G'{t) in the FPT problem for a Brownian motion. Writing 

^'(z) := -^exp {^—\^, the derivative g satisfies a Volterra integral equation 

of the first kind of the form 

This and other such integral equations can be used to find g numerically. 

A. Shiryaev is generally credited with introducing the IFPT problem in 1976 
(we have not been able to find an explicit reference). Most authors have inves- 
tigated numerical methods for finding the boundary. Details can be found in 
[HWOOi IHWOli IIK021 [ZS09] . It is shown in jAZOl) that for sufhciently smooth 
boundaries the density u{x, t) and the boundary b{t) are a solution of the following 
free boundary problem 



(1.8) 



ut{x,t) 
u{x, t) 
u{x, 0) 

G{t) 



^(o-^u)^;^ - {tJ'Ujx, X > b{t), i > 0, 

0, a; < b{t), t > 0, 
f{x), xeR, 



u{x, t) dx, t > 0. 



where / is again the probability density of Yq. The existence and uniqueness of a 
viscosity solution of (|1.8|) is established in [ CCCSll] along with upper and lower 
bounds on the asymptotic behavior of b. This paper also shows that this b does in 
fact produce a boundary that gives the survival function G. To our knowledge it 
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has not be proven that a strong solution to the system (jl.Sp exists, nor that there 
is a smooth b solving the IFPT. 

A variation of the IFPT is studied in IDPllj . There the barrier is fixed at zero 
(i.e. 6 = 0) and it is the volatility parameter cr(-, •) that is allowed to vary. The 
authors show that this problem admits an explicit solution for every differentiable 
survival function. 

2. Global Existence and Uniqueness 

Suppose for the remainder of this paper that Yj :— Yq + Bt where {Bt)t>o is a 
standard Brownian motion and Yq is a random variable, independent of B and with 
density /. In this case, (|1.2p is 

ft 



Git) 



E 



exp 



A / tp{x + B, ~b(z))dz 



fix) dx 



which, by time reversal, becomes 



Set 
(2.1) 



G{t) 



u{x, t) 



exp - / ^(x + Bt-, - b{z))dz fix + Bt) 



dx. 



E 



exp 



^ix + Bt-, - biz))du] fix + Bt) 



That is, u is the sub-probability density of Y killed at the random time r. It is well 
known that the function u is the unique solution of the PDE 

Utix, t) — -u^xix, t) — Xipix — bit))uix, t), X & R, t > 0, 
uix,0)^fix), xeR. 
Any solution to this PDE satisfies 



(2.2) 



lim uix,t) — lim Uxix,t)=0, t > 0. 

a:— ^±oo a:— ^±cxD 



Our question as to whether we can find a "barrier" b giving us the survival function 
G is now equivalent to whether the system 

Utix, t) ~ -Uxxix, t) — \'4'ix ~ bit))uix, t), X £R, t > 0, 



(2.3) 



uix,0) = fix), x&R, 
uix,t)dx = Git), t > 0, 



has solutions (m,6). Differentiating the third equation from (|2.3I) with respect to t 
and then using the first equation together with an integration by parts, we get that 

(2.4) -git) = \ / ipix-bit))uix,t)dx, 

where we recall that git) — G'(t). A second differentiation in t followed by another 
integration by parts yields 

g'it)-^^ / 'ip'^ix-bit))uix,t)dx 



A / ^xix -bit))uix,t)b'it)dx 
— A/2 / Tpix — bit))uxxix,t) dx 
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(2.5) 



A / il)x{x-b{t))u{x,t)b'{t)dx 
Jr 

+ A/2 / Tp^{x — b{t))ux{x,t) dx 
Jr 

A / il>x{x -b{t))u{x,t)b'{t)dx 
Jr 

- A/2 / tp^.^ix ~ b{t))u{x,t) dx. 



Note that (|2.5p may be rearranged to give an ODE for b of the form b'{t) — 
Q{b{t),t)^ where the function Q is constructed from the function u (which, of 
course, depends in turn on b). Re- writing this integral equation in the form 
b{t) = b{0) + Q{b{s), s) ds leads to the following theorem, our main result. 

Theorem 2.1. Suppose the following. 

• The survival function G is twice continuously differentiable with first and 
second derivatives g and g' and < —g{t) < XG(t) for all t > for some 
constant A > 0. 

• The initial density f satisfies Jg^f{x)dx = 1, f{x) > for all a; G 
/ G C^(R), and the functions /, /', /" are bounded. 

• The function ip is non-increasing and belongs to C'^(R), and for some h > 0, 
^{x) = 1 for x < —h and ip{x) — for x > h. 

Then, there exists a unique continuously differentiable function b such that the 
following three equations hold 



(2.6) 
(2.7) 



G{t) 



[ E 


exp(-A^ 


Jr 





A / ip{x + Bu — b{u))du 



lit) = A 



fix) dx, 



exp -A / i/j{x + Bu- b{u))du ^:{x + Bt ~ b{t)) 



fix) dx, 





■^2(3. j^B,- b{s))e-^ /o ^(x+B^-b(r))dr 


fix) dx 




^xi 


X + Bs- b{s))e-^ i,{x+B.~b{r))dr 


fix) dx 



+ 



ds 



and 

bit) = 6(0) 

(2.8) 
for all t > 0. 

Proof. From now on we assume for ease of notation that A = 1. The modifications 
necessary for general A are straightforward. The proof will be via a sequence of 
lemmas, all of them assuming the hypotheses of Theorem 12.11 (with A = 1). We 
start with the following simple observation. 





^xxix + Bs- h{s))e-^ K ^(^+B^-b(r))dr 


fix) dx 




^P^{X + Bs- b{s))e-^ ^0 ■^^'■'+Br~b{r))dr 


fix) dx 



Lemma 2.2. Suppose that 



Git) 



u[x, t) dx 
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for some continuous function u : M.X — > M such that u{x, t) > for a; G M, t > 0. 
Then, for each t > there exists a unique b{t) € R such that 



Proof Set 



Then, 



-git) 



F{t,z) 



■ip{x — b{t))u{x, t) dx. 



^{x — z)u{x, t) dx. 



lim F{t,z)^ / u{x,t)dx = G{t), 



and, by assumption, 



lim F{t,z)^0, 

2— )- + 00 



< -g{t) < Git). 



Furthermore, F is continuous and strictly decreasing in z. So, by the intermediate 
value property, we can find a unique b{t) £ M such that F{t, b{t)) — —g{t). □ 

Lemma 2.3 (Global Uniqueness). Suppose there exist continuous functions 6i,&2 
such that equations (|2.6p . (j2.7p and (|2.8p are satisfied for b ~ bi and b ^ b2. Then, 
bi{t) = 62 (i) /or allt>0. 

Proof. Recall that we are assuming A = 1 to simplify notation. 

Suppose that bi and 62 are two continuous solutions of ()2.6|) . ()2.7|) and (j2.8|) . It 
follows from Lemma O and dSj]) that 61 (0) = 62(0). Set V inf{i > : bi{t) ^ 
62 (i)} and suppose that V < 00. 

Define / : R ^ M by 



f{v)dy 



E 



+ ^ M, i = 1,2, by 



where 6(t) = fei(t) — ^2(i) for < t < T^. Define functions b. 
hit) ^ bi{V + t),t> 0. Then, 5i(0) = 62(0) = 6(F), and 

6,(i) = 6,(0) 





^2(^ + B« - 6^(s))e- ^0 dr 


f{x) dx 




?/'x(a; + Bs - 6,(s))e-/o V'(^+B.-b,(r))dr 


f{x) dx 



1/2 /rE 


^PrAx + Bs - b,{s))e-io■^'^-+B.-i^{r))dr 


f{x) dx 


/rE 


i:^{x + B, - biis))e- fo ^^^+Br-b^(r))dr 


fix) dx 



ds. 



Fix e > and set 



K := min inf 

1=1,2 0<s<e 



Mx + Bs - 6,;(s))e--''o fix)dx > 0. 



By the triangle inequality, for < t < e, 

\bi{t) - b2{t)\ < I + 11 + III, 
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where 



I:=K-' f\g\8 + V)\ [ E 
Jo Jr 



f{x) dx ds, 



II := K-^ 

X / E 



E 



+ B,- b2{s))e- 'A(^+B.-fe(r))dr 



i)^{x + B,- h{s))e- ^0 M^+Br-h{r)) dr 



f{x) dx ds 



+ K- 



11 

Jo Jr 



E 



^P^x + B, - bi{s))e- 



- tP(x+Br-bi(r))dr 



.^■^i^X + Bs- b2{8))e- ^i^+Br-Ur)) dr 



Jr l 



Rx) dx 



x-\- Bs- bi{s))e 



■J^^{x+B^-bi(r))dr 



f{x) dxds, 



and 



III := -K-^ 
2 



X / E 



E 



iIj^^{x + B, - 6i(s))e-/o V'(x+Sr-6i(r))dr 1 f(^x)dx 



I 

Jr 



Xl)^{x + B,- b2{s))e- fo Hx+B^-h{r)) dr 



- + - h (S))e- -^o V'(x+S,-6i(r)) dr 



f{x) dxds 



+ 



)-K-^ f I e[ Vx.(a; + B«-6i(s))e--/'o'^(-+Sr-6iW)dr 



^PUX + b2{s))e- /o ^iM-b.ir)) dr 



f{x) dx 



X / E 



V(x+Sr-6i(r))dr 



f{x) dx ds. 



ipx{x + Bs- 6i(s)) 
Consider the integrand in /. Note that 

i^^ix + B,- b2{s))e- ^0 ^i^+Br~Ur))dr b^(s))e' ^0 ^i^+Br-bi(r))dr 

+ e-/o + B,- b2{s)) - + - h{s))\ 

< ||V'x||i,~s||V'x||i,<~ sup \b2{r) - bi{r)\ 

0<r<s 

+ ||V'xx||l~ sup \b2{r) - bi{r)\. 

0<r<s 

Similar arguments for the integrands in // and /// using the boundedness and 
global Lipschitz properties of ip, ipx, and V'xx establish that, for a suitable constant 
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0<s<t 



sup l&i 62(5)1 <C / sup \bi{r)~b2ir)\ds 



0<r<s 



for < i < e. It follows from Gronwall's inequality that bi{t) = b2{t) for < i < e, 
and so bi{t) ~ b2{t) for < t < V + e, contrary to the definition of V and the 
assumption that V is finite. □ 

Lemma 2.4 (Global Existence). Define S to be the supremum of the set ofT such 
that the equations (I2.6p . (|2.7p and (|2.8p have a continuous solution on [0,r]. Then, 
S = +00. 

Proof. Suppose to the contrary that S < +00. From Lemma 12.31 the equations 
have a unique solution on [0, S). By time- reversal, equation (|2.6|) is equivalent to 



(2.9) G{t) 



E 



exp 



ij{x + Bt-u ~ b{u))du fix + Bt) 



dx. 



Similarly, p.7|) is equivalent to 

-9{t) 

(2.10) 



exp ( - / ^j{x + Bt-u - b{u)) du ) i}{x - b{t))f{x + Bt) 
¥ovQ<t<S put 



dx. 



(2.11) 



u(a;, t) := E 



exp - / iix + Bt-u - b{u))du fix + Bt) 



Consider ti < t2 < . ■ . t S . It follows from the continuity of the sample paths of 
B that as i„ f S 



exp 



^j{x + Bt^-u - Ku)) du fix + BtJ 



exp 



^ix + Bs-u - biu)) du fix + Bs) 







almost surely for each x € M, and so 

rS 



uix, tn) E 

Because 



exp 



^Pix + Bs-u - biu)) du fix + Bs) 



uix, S). 



uix,t)<E[fix + Bt)], 
it follows from dominated convergence that 

uix,S)dx~\im / it(a;, i„) da: = lim G(i„) = G(S'). 



Also, 

lim / ipix — bit,i))uix, tn) dx — ~ lim (;(f„) — —giS). 

Because < -giS) < GiS) and 

uix,S)>e^'^E[fix + Bs)]>0, x G R, 
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there is, by Lemma [221 a miique 6* G K. such that 

ip{x — b*)u{x, S) dx — —g{t). 



We claim that b{tn) b* . If this was not the case, then, by passing to a 
subsequence we would have converging to some other extended real c and 

hence, by dominated convergence, 

~9{t) = -limg(i„) 



lim J ^p{x — b{tn))u{x,tn) dx 
ip{x — c)u{x, S) dx, 

contradicting the definition of b* (where we used the natural definitions ip{—oo) :~ 
l,V'(+oo) :— 0). Using dominated convergence in (|2.8p we get that there exists a 
continuous b such that all three equations hold on [0,S]. 

All we need to do now is show that we can extend the existence from [0, S] to 
[0, 5 + (5] for some ^ > 0. This amounts to proving existence on [0, 6] starting at 
a different initial condition - replacing the original probability density / by the 
density of the probability measure 



E 



exp 



ip{x + Bu- b{u)) du] ,Bs ^ 



fix) dx/G{S). 



This will follow if we can establish the local existence, that is the existence for some 
(5 > 0, of a solution of the following PDE/ODE system 

1 . 



ut{x, t) — -Uxx{x, t) — 'ip{x — b{t))u{x, i), X G M, < t < (5, 
u{x, 0) 
6(0) 



b'{t) 



u{x,S)/G{S), X eR, 
6(5), 

jgjS + t)+ g'jS + t))/GiS) - /r[V''(x - bjt)) - ^(x - b{t))\u(x, t) dx 
Jk''Px{x - b{t))u{x,t) dx 
1/2 Jg^ipx{x - b{t))ux{x,t) dx 



Jg^ipx{x - b{t))u{x,t) dx 



-,0<t<6. 



We note that the expression for b'{t) is not the analogue of the one for b'{t) that 
arises immediately from differentiating (|2.8p . which in turn arose from rearranging 
(j2.5|) and integrating. However, adding = J^^j{x — b{t))u{x,t) dx — g{t) to the 
right-hand side of (I2.5P and then solving for b'(t) leads to an expression of this 
form. Note that 

iix + Bs-u~b{u))dii\ f{x + Bs) 



u{x, S) 



exp 



>0, 



and, by dominated convergence, that u{-,S) G C^(K) with ||m(-, S')||ioc(]g'), 
l|wa;(-, S')||j;,oo(j{), ||uxa;(-, S')llioo(R) aU finite. Therefore, we can apply Theorem |3JJ| 
below to get that there is a time 5 > and a unique pair u, 6 satisfying the 
PDE/ODE system above with u twice continuously differentiable in a; on R and 
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once continuously differentiable in t on [0,6], i.e. u G C'^{M.)Cf{[0,6]), and with 
b € (^^([0, S]). Thus, we have proven that we have a unique continuous b satisfying 
equations (|2.6p . ()2.7|) and (|2.8|) on [0, S* + (5]. This contradicts the maximahty of S. 
As a resuh, S = oo and we are done. □ 



This completes the proof of Theorem 12. II 



□ 



As a corollary we get the global existence and uniqueness of the PDE/ODE 
system. 



Corollary 2.5. Suppose that the conditions of Theorem \2. l\ hold. Then, the system 

ut{x, t) = ^Uxx{x, t) - i}{x - b{t))u{x, t), 
u{x,0)^f{x), x£R, 

-5(0)= / ^{x-b{0))f{x)dx, 

JR 

g{t) + g'{t) - - b{t)) - ij{x - b{t))]u{x, t) dx 



(2.12) 



b\t) = 



/k V-'a: {x - b{t))u{x, t) dx 

1/2 J^^pxix - b{t))ux{x,t)dx 



iK^xix - b{t))u{x,t) dx 
has a unique solution {u,b) e Cl(R)C}{R+) x C}{R+). 



t > 0, 



3. Local Existence and Uniqueness 



We now consider the PDE/ODE system (|2.12l) . We have already used the stan- 
dard notation Fx and Fxx to denote the first and second derivatives of a function 
F of one variable or the first and second partial derivatives with respect to the 
variable a; of a function F of several variables. Because we repeatedly deal with the 
function {x,t) i— >■ ip(x — b{t)), it will be convenient to recycle notation and define a 
function ipb -.Rx M+ — > R by ipbix, t) — ^{x — b(t))}. We will then set 'ipx,b '■= dxtpb 
and '>pxx,b ■= dxx^b- We will continue to use the notation ipx and ipxx with its old 
meaning, but there should be no confusion between the different objects ipi, and tpx- 
Similarly, we set (j) ■= ij;'^ — ip — — ^) and put (j)i,{x,t) = (f){x — b{t)). Lastly, 

for two functions f,g define {f,g) = j^f{x,t)g(x,t)dx. 

In the notation we have introduced, we wish to consider the system 



(3.1) 



"^tix, t) = ]^Uxx{x, t) - ip{x - b{t))u{x, t), X e 
u{x,0)^f{x), xeR, 
&(0) - bo, 

g{t)+g'{t) - {(f>b,u) - l/2{il;x,b,Ux) 



b'{t) 



l,t>Q, 



t > 0, 



{lpx,b,u) 

for some bo G R. (In the proof of Theorem 12.11 we choose bo to satisfy —g{Q) = 
/gt/'(a; — bo)f{x)dx, but we may take an arbitrary value for bo and still obtain a 
local existence and uniqueness result.) 

We have assumed in the statement of Theorem 12.11 that f E (K) and tp G 
C^(R) with UWl^ = 1, IIV'IIl- =: B, Uxx\\l^ C, and Wi^xxxh^ =: F for finite 
constants B,C,F. Furthermore, we have assumed for some h > that ^Pi^) — 1 
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for X < —h, that tlj{x) — for x > h, and that ip >Q and ipx < for all a; G M. Set 
/r d-x =: D and note that < I? < oo. It is immediate that ||0||/^oo < 1 and 

= " 2V')||l°° < \\fpx\\L°° = B. Moreover, the functions and (l)x are 

supported on [—h, h] and < \<t>(x)\ dx =: E < oo. 

Definition 3.1. For T > 0, let {C'^ , \\ ■ ||t) be the Banach space consisting of pairs 
of functions {u,b) such that u G C^{R)Cti[0,T]), b e C{[0,T]) and 

||(u, 6)||t := ||w||l~(r)l~([o,t]) 

j-g 2) + ll"^IU~(R)Lr([0:T]) + II"^:i=IIl~(R)L-([0,T]) 

+ II^IIl-([o,t]) 

< oo. 

Definition 3.2. Given constants M,N,P,A,L > 0, foo e M and T > 0, define the 

closed subset ^MNPALbo 

II"IIl~l~([o,t]) < M, 

ll"a:||Lr([0,T])L~ < N, 
(3.3) \\Uxx\\L^i[0,T])L^ < P, 

b{0) = 6o, 

I16||l~([o,t]) < A/2, 

inf t) > L 

xe[-A,A],te[o,T] 

The following is the main result of this section. 

Theorem 3.3. Suppose that the assumptions of Theorem \ 2.1\ hold. Suppose also 
that the constants M, N, P, A, L > and 6o G M are such that 

• \bo\ < A/4, 

• fix) >4:L>0forxe [-A, A], 

• ||/||l="(r) <M/2, 

• II/x||l-(r) < N/2, 

• II/3:2:||l=(R) < -P/2. 

Then, for T > sufficiently small there is a contractive map $ : ^'MNPALbo ^ 
^MNPALbo defined by ^{v,b) = {u,c), where 

ut{x, t) = —Uxx{x, t) — ^{x — b{t))v{x, t), X £ M.,t > 0, 



u{x,0)^f{x), xeR, 

f,,. gjt) + g'jt) ~ {<f>b,v) - l/2{i;x,b,vx) r, ^ , ^ rj. 
{Vx,b,v} 

c(0) = bo. 



(3.4) 



We will prove Theorem 13. 31 in a series of lemmas. Each lemma will assume the 
hypotheses of Theorem 13.31 and the bounds established in the previous lemmas. 



Remark 3.4. Since / is continuous and positive, for any A > there exists L > 
such that f{x) > 4L for x G [—A, A]. Therefore, we are not restricting the possible 
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values of 6(0) by the above assumptions. We will also assume without loss of 
generality that h < A/ A. 

Lemma 3.5 (Boundedness of u). Suppose that {u,c) — with {v,b) G 

^TiNPALbo ■ Then, there exists a time T > such that 

WuWl^l^Ho.t]) < M. 
Proof. Using Duhamel's formula (see ()8.2p ). 



\u{x,t)\ 



G{y,t)J{x-y)dy- 



G{x -y,t- s)i!c(s) iy)viy, s) dy ds 



< 



/ G{y,t)f{x-y)dy+ / / G(x - y,i - s)|?/'c(s)(y)||w(2/, s)| 

Js. Jo JR 



ds 



< M/2 / G{y,t)dy + M / / G{x - y,t - s) dy ds 

< M/2 + Mt 

< M 



when t < h, where 



G'(x,t) := ^=e" — , xeR,t>0. 



□ 



Lemma 3.6 (Boundedness of Ux)- Suppose that {u,c) — ^{{v,b)) with {v,b) £ 

a time T > such th 

II^2:||l~([0,T])L~ < ^■ 



^TiNPALba ■ Then, there exists a time T > such that 



Proof. Since Ux solves 



dr. 



dt - 1 Ux 



-l/jx.cV - IpcVx, X & R, t > 0, 



Uxix,0) = fxix), 

we have via Duhamel's formula that 



\Ux{x,t)\ 



G{y,t)fx{x - y)dy 



< 



+ / G{x- y,t- s){-ipx,cV- il'cVx)iy,s)dyds 



< / G{y,t)\fx{x - y)\dy + / / G{x - y,t - s)\'4>x,c\\'"iy, s)\ dy ds 



G{x -y,t- s)\il;{y - c{s))\\vx{y, s))\ dy ds 



N 



+ MB 



G{x — y,t — s) dy ds + N / / G{x — y,t — s) dy ds 



N 



< — + MBt + Nt. 
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Thus, 



N 



MLr{[o,T])Lr < Y + (MB + N)T<N 
whenever T <T'^^ where 



TV 



□ 



2 (MB + TV) 

Lemma 3.7 (Boundedness of Ux^)- Suppose that {u,c) — ^{{v,b)) with {v,b) G 

a time T > such th 

Uxx\\l-^{[0,T])L^ < P- 



^A/TVPALfco ■ ^^s*^; there exists a time T > such that 



Proof. Noting that Uxx solves 
Ot 



•XX = -i'xx.cV - 2lpx,cVx - tpcVxx, X E R, t > 0, 

analogous manipulations to those from Lemma 13.61 yield the result. 



□ 

Lemma 3.8 (Lower bound for u and boundedness of c' and c). Suppose that 
{u,c) — $((«,&)) with {v,b) G ^'A/[NPALbo' Then, there exists a time T > such 
that 

(3.5) u>Lonxe[-A,A],te[0,T], 
and c{t) e [-A/2, A/2] for t e [0, T] . 

Proof Recall that 6(0) £ [-^4/4, A/4]. Then, it is immediate that 



(3.6) 



ipxix — b{t))v{x, t) dx 



'4'x{y)v{y + b{t))dy 



>DL, te[0,T], 



because on the support [~h,h] of tpx we have y e [—h,h] C [—A/A, A/A] which 
together with the bound on 6(t) implies t/+6(i) G [—A, A]. Therefore, v{y+b{t)) > L 
for i e [0,T] which, since < 0, yields 



MyHy + b{t)) dy<L j My) dy = -ld < o, te [o, r]. 

unds thi 
c'{t)\ < 



We see from these bounds that 

snp[a4\9 + 9'\)+ME + ND/2 



LD 



and, by integrating, 



suprn ti(|5 + g'\) + ME + ND/2 
\cit)\ < |c(0)| + '-t. 

Thus, there is T > such that for t e [0, T], 

\c{t)\e[-A/2,A/2]. 
Using the assumptions, equation (|8.2p gives 



u{x, t) 



G{y,t)f{x-y)dy 

■x+A 



G{x -y,t- s)V'c(s) {y)v{y, s) dy ds 



> AL 



G{y,t)dy~M / / G{x - y,t - s) dy ds 
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px+A 



> AL I G{y,t)dy-Mt. 

Jx-A 



Ix—A 

li < X < A then a; — ^ < and x + A>A>Qso for small enough t we have 



G{y,t)dy> I G{y,t)dy> 

-A 



1 



If —A < X < then x + A>0 and x — A < —A < 0. So, for small enough t, 

rx+A 



f G{y,t)dy> f G{y,t)dy>\. 

Jx-A J-A •J 



Therefore, there exists a time T > such that whenever t e [0, T] and x G [-^A, A], 

4 

u{x,t) > -L-Mt 

O 

> L. 

□ 

Lemma 3.9. For a sufficiently small time T > 0, the set TMiVPALbo mapped 
into itself by $. 

Proof. The above lemmas provided the necessary bounds. Now, note that if we 
start with {v,b) € ^'AjNPALbo' then we first get the function c from the last two 
equations in (j3.4l) by simply integrating. The integration is well-defined because 
the denominator is bounded in absolute value below by DL > and the numer- 
ator is bounded above. Thus, c G C^([0,t]). Next, having c in hand we get the 
function u from the first two equations of ()3.4p . We note that, by Duhamel's 
formula, the function u has actually more than the desired smoothness, namely, 

7. GC2(R)ci ([o,r]). □ 

Lemma 3.10. Suppose that (ui, 6i), (w2, ^2) G ^JjNPALbo- {ui,ci) = 

61)) and (u2,C2) — <i>((f2, &2))- For any e > there exists T > such 

that 

(3.7) \\c2 - ci\\l^(^io,t]) < e|l(w2,&2) - (ui,&i)||t- 

Proof. Note that the functions ci , C2 satisfy 



(3.8) 



/ . ,N 9{t) + g jt) - {4>b, , vi) - i^ii^xM ^ dxvi) , ^ „ 

Ci[t) = py- ^ , t>0, 

C2(t) = j^y- ^ , t>0. 



Subtracting the two equations gives 

4(0 -cUt) = [g{t)+m ( ^t-"-"^^ - i^--"-",^^ + ^t7^-"^^ " i^^M'-^) ) 

{{(l^bi , Vl) - {(t)b., , Vi)) {ljja:,b2 , V2) , ((062 > «l) ~ ('^62 ) ^2)) (V'a:,&2 > ""2) 



+ 



(^x,bi , Vl){lpa:,b2 , V2) {ipxM > «l) (^a:,62 ) ^2) 

(('/>b2 7 ^^2) - {(|}b^ , ^2)) ('/>b2 ^ ""2) ^ ((0bi , V2) - {<Pbi , Vl)) (0b2 . V2) 

(V'x,bi , Vl){i>x,b2 , ^2) (V'x,bi , Vl){lf>x,b2 , ^2) 

{{■^Jx,bi,dxVi) - {ipxM^^xVi)) {l}} 

x,b2 ; '^2/ 



16 BORIS ETTINGER, STEVEN N. EVANS, AND ALEXANDRU HENING 

+ 
+ 
+ 



((V'X,62. 


dxVi) 








i;,fcl,^^l)('0ri:,b2:'^2) 


((V'x,62, 


V2) - 


(V'i;,6i,«2)) {llJx,b2^9xV2) 






6i,Wl>(V'a;,h2.«2) 




V2} - 


{ipxM'I'l)) {lpx,b2,dxV2} 



2{tpxM''"i){'^x,b2^'"2) 

Using the fact that the functions tp, ipx and (f) arc Lipschitz, that Vi and V2 are 
bounded, and that their first derivatives are bounded, we find that 

„/ /„ / SUP[o,T] Ifl + ff'llkl -^^2||l~L~([0,T]) 
I|C2-Ci||l~([0,T]) < 

suP[o,T] \9 + 9'\MC{A + 2h)\\b2 - 6i||l~([o,t]) 



+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 



DM^B{A + 2h)\\b2-b4Lrilo,T]) 

DME\\V2 - WllU^Lf ([0,T]) 

EM^B{A + 2h)\\b2-b4Lrao,T]) 

ME'^\\V2 — Vl\\ {[0,T]) 

NMDCjA + 2h)\\b2 - 6i||Lr([o,r]) 

MD'^\\dxV2 - 9a;Wl||L~I,oo([0,T]) 

NMDCjA + 2h)\\b2-bih^^^o,T]) 
2L2D2 

ND^\\V2 - fl||L~L~([0,T]) 



2X21)2 

Integrating and recalling that ci(0) = C2(0) = 60 leads to 



= |C2(t)-Ci(t)-(C2(0)-Ci(0))| 



/ {C'2{s)-C[{s))ds 

Jo 

< I \C'2{S)-C'^{s)\ds 

Jo 

< ^l|c2-ci||L~([0,t])- 

Hence, 

\\C2 - Cl|iL~([o,T]) < ^114 - c'i|lioo([o -rj), 

and by the above bound on ||c2 — c'i||/,cc([o,t]) for any e > we can choose T small 
enough that 

||C2 - ci||l«>([o,t]) < e\\{v2,b2) - {vi,bi)\\T- 

□ 

Lemma 3.11. Suppose that {vi,bi), {v2,b2) S ^'uNPALbo- {ui,ci) = 

^{{v\,b\)) and (u2,C2) = $((v2,62))- For any e > there exists T > such 
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that 

(3-9) \\u2 - uilU^irdO'^]) - <^ll("2,&2) - (wi,6i)|lr- 

Proof. The following equations hold 



(3.10) 



dt-^jui^-xl:{x-ci{t))vi, xgM, t>0, 

dt - U2 = -1p{x - C2{t))v2, X & R,t>0, 

ui{x,0) = f{x), xeR, 
U2ix,0) = f{x), xeR. 

By Duhamel's formula we have 

(3.11) ui^G^ifSt^o) + G*i-^p,,vi) 
and 

(3.12) u2^G*{fSt^o) + G^i-^c,V2), 

where we recall that * denotes convolution on x R. Subtracting the two equations 
gives 

Ui~U2^G* {{tpc2 - ^ci )V1 + ^C2 {^2 -vi)). 

Bounding in terms of the sup norm and using the fact that 

\i:{x - Ci{t)) - ^P{X - C2m < ||^A.|||ci(t)-C2(i)|, 

we have 



\ui{x,t) - U2ix,t)\ < 



G{x-y,t~s)\Tpc^{y,s) - ipc2{y, s)\\vi{y, s)\dy ds 



G{x -y,t~ s)\'ipc2iy, s)\\v2{y, s) - vi{y, s)\ dyds 







< IIV'2;||L~||t^l||L~L~([0,T])||ci - C2\\L^t 

+ llV'Ik-ll^^l - '"2||L-L-([0,t])i 
= BM\\ci - C2\\L^t + \\vi - W2|lL~L-([o,t])i- 



Thus, 



- U2\\l^L^{[0.T\) < B\\ci - C2|1l~T + 11^1 - V2\\l^L^{[0,T\)T. 

so for small enough T we see that p.9p holds. □ 

Lemma 3.12. Suppose that (wi, 6i), (f2, ^2) €E ^\,iNPALbo- (ui,ci) = 

<i>((tii, 61)) and (u2,C2) = *&((f2,fe2))- For any e > there exists T > such 
that 



(3.13) 



\dxUi - dxU2\\L^L^([o^T]) < e|l(i'2,&2) - (wi,6i)|1t- 
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Proof. Differentiating p.lO|) with respect to 



(3.14) 



x£R,t>0, 



9x 



dt —j dxU2{x,t) = -i^x,C2{x,t)v2ix,t) - li)c2{x,t)dxV2{x,t), 

X e K,t > 0, 

d^ui{x,Q) ^ f^{x), .tgM, 
d^U2{x,Q) ^ f^{x), .tgM, 
Via Duhamel's formula, 

(3.15) d^ui = G * (/,(5t=o) + G * - ci{-))vi - ^(- - C2{-))d^vi) 
and 

(3.16) d^ui = G * ifJt=o) + G * i-M- - C2{-))V2 - - c2{-))d^V2). 
Subtracting and rearranging, 

{dxUi- dxU2){x,t) ^ j j G{x-y,t- s)[iljx^c2V2iy,s)-^x,ciVi{y,s)]dyds 

+ / G{x-y,t- s)[ilJc2dxV2{y,s)-ipcidxVi{y,s)]dyds 
Jo Js. 

G{x ~y,t~ s)['iljx^c2V2{y,s) - V'a;x2«i(y, s)] dy ds 
G{x -y,t- s)[tp^^c2Vi{y,s) - tpx,ciVi{y,s)] dyds 
+ 11 G{x ~ y,t ~ s)[ipc2dxV2iy,s) - ipc2dxVi{y,s)]dyds 



G{x -y,t- s)[%l:c2dxVi{y,s) - iic^dxVi{y,s)] dyds. 



Using estimates similar to those in the proof of Lemma 13.111 

\\dxUi- dxU2\\L-^Lf=-([0,T]) < -BM||U2 - Ul||L~Lj-([0,Tl)r 

+ GM\\c2 - ci\\l^{1(),t])T 

+ \\dxV2 - (9a;Wl||L~L~([0,T])r 

+ BN\\c2 - ci||l~([o,t])7' 
= BMT\\v2- vi\\l^l^(^io^t]) 

+ {GM + BN)T\\c2-ci\\Lrao,T]) 

+ T\\d:^V2 - 9a;Wl||L-L~([0,T])- 

SO for T small we recover p.l3p . □ 

Lemma 3.13. Suppose that (wi, &i), (w2, ^2) € ^fiNPALbo' (wi,ci) — 

*i'((t^ii ^1)) o^nd {u2,C2) = '^{{v2ii>2)). For any e > there exists T > such 
that 



(3.17) 



\\dxxUi - 5j;3;U2||l^l~([o,t]) < e||(w2,fe2) - (ui,6i)|1t- 
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Proof. Differentiating p.lOp twice with respect to x 
3 



(3.18) 



dt 



3 
~2 



dxxUl = -■)pxx,ciVl - 2'>px,cidxVi - IpcidxxVl, 

xeR,t>0, 
X eR,t>0, 



3xxUi{x,0) ^ fxx{x), a; e R, 
3xxU2ix,0) ^ fxx{x), a; e R. 

Duhamel's formula and similar manipulations to Lemmas 13.111 and 13.121 give 

\\dxxUi - dxxU2\\L^L^{lO,T]} < C'A^ II ^^2 " Wl || L~ ([0,T])L~?' 

+ FM\\C2 - Ci||i=.o([o,T])7' 

+ 2B\\3xV2 - 3xVi\\L^L^(^[a^T])T 
+ 2CN\\c2 - ci||l~([o,t])T 

+ \\dxxV2 - dxxVl\\L^L'i°(lo,T])T 

+ BP\\c2 - ci||io=([o y])T 
= CMT\\v2- vi\\l^l^(^[q^t]) 

+ 2BT\\dxV2 - dxVi\\L^L^(^in,T]) 

+ T\\3xxV2 - dxxVl\\L^L^([o^T]) 

+ {FM + 2CN + BP)T\\c2 - ci||l~([o,t])- 
so when T > is small (PTT)) holds. □ 

Theorem 3.14. [Local existence and uniqueness] Suppose that the conditions of 
Theorem \2.1\ hold. Then, there exists a time T > such that the system 



ut{x, t) = -Uxx{x, t) — '\fi{x — b{t))u{x, t), X &R, t > 0, 



u{x,0) = f{x), xeR, 

g{t) + g'it) ~ {(pb, u) ~ l/2(^^_b, Ux) 



b'{t) 



{lpx,b,u) 



t > 0, 



6(0) = 6o, 

has a unique solution {u,b) £ Cl{R)C^ {[0,T]) x C^{[0,T]). 

Proof. Note there exist strictly positive constants A, M, N and P such that bo e 
h4'4]' /(^) >L>0, when x € ||/|Uo.(r) < Af, ||/,||loc(k) < N/2, and 

||/i^||l=(r) < P/2. Putting all the estimates from the above lemmas together we 
have that, if < e < 1 is fixed, then for T > small enough 



U2,C2) - (ui,ci)|| < e||(w2,62) - {vi,bi)\ 

T : 

a contraction. Since FT, 



Thus, there exists a T > such that the map <i> : 



MNPALbo ~^ ^ MNPALbo 

MNPALb ^ closed subset of the Banach space £^ , the 
Contraction Mapping Theorem gives that there exists a unique fixed point, that is. 
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a pair {u,b) e C^{R)Ct{[0,T]) x C([0,T]) with 6(0) = bo such that 



(3.19) 



ut{x, t) = -Uxx{x, t) - V'(x - b{t))u{x, t) 

git) + g'{t) ~ {cbb,u) ~ l/2(i^,^b, u,) 



b'it) = 



6(0) = bo. 

We can now argue that our fixed point (u, b) has more smoothness than it seems 
a priori. The third equation in (I3.19P imphes that b must be continuously difFeren- 
tiable with a bounded derivative. This, together with the first equation from (I3.19P 
then tells us that u has a continuous derivative in time. Therefore, we must have 
iu,b) € C^^iR)C}{[0,T]) X C\[0,T]). □ 



Corollary 3.15. Assume the hypotheses of Theorem \3. 14\ and the extra conditions 

G(0)= / f{x)dx, 

JR 

(3.20) 



-5(0)= / ^P{x-biO))fix)dx, 
Jr 

0<-g{t)<G{t), te[0,T]. 
Then, there exists a time T > such that the system 

Ut{x, t) = ^u^xix, t) - ^{x - b{t))u{x, t), X el 
u{x,0) = f{x), xeR, 

G{t) = / u{x,t)dx, t e [0,T], 



0<t<T, 



has a unique solution {u,b) : M x [0,r] 
andbeC\[0,T]). 



Furthermore, u € Cl{R)C^ {[0,T] 



Proof. First note that by Lemma 12.21 we have that 6(0) is uniquely determined. 
From Theorem 13. 141 we have that there exist unique u, 6 with u £ C^(M)C(^([0, T]) 
and 6 e C^([0,T]) satisfying the PDE and having everywhere in [0,r] 

^ git)+g'it) - {^b,u) - i/2{iPx,b,ux) 

Set F{t) :— G{t) ~ J^u{x,t) dx and note that the first two conditions from 
(P:^ yield, together with the PDE, Ft{0) = F{0) = 0. The function F belongs 
to C^([0,r]) and Ft belongs to C([0,T]). The above equation for 6' is equivalent, 
after using the PDE, to 

Fu{t) - Ft{t) ^ 0, te[0,T]. 

Integrating and using the fundamental theorem of calculus, we get 

Ft{t) - Fit) = FtiO) - F(0) = 0, te [0,T]. 

The unique solution to this differential equation is F(t) — Ce* for some constant 
C e M. This together with F(0) = yields F{t) = for f e [0, T]. Thus, 



Git) 



uix,t)dx, t e [o,r]. 
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Then, taking a derivative and using the PDE, 

-g{t)= I i}}{x -h{t))u{x,t)dx, tG[0,T]. 
Because |'0(a;)| < 1 for a; e M, "i/; = for x > h and u{x, t) > we see that 
0< / i:{x - b{t))u{x,t)dx = -g{t) < / u{x,t) dx = G{t). 



□ 

4. Discontinuous killing 

Next, we consider the existence of a barrier when kiUing is done non-smoothly. 
That is, we ask whether there exists a function b such that, for a given G 



(4.1) G{t) / E 



cxp ^- J l(_oo,o](a; + Bu- b{u)) diij f{x) 



dx 



l^{x)= / ^((y-e)/e))(l/6)d2/ 

J X 

^,{x) = / 0(y/e))(l/e)dy, 



Note that l(_oo,o](2^ + ~ ^(w)) du is the time during the interval [0,t] spent 
by a Brownian motion started at x below the barrier b. 

Theorem 4.1. There exists a function b such that, for a given, twice continuously 
differentiable G satisfying < —g(t)/G{t) < 1, t > 0, equation (|4.ip holds for all 
t > 0. 

Proof. Let be a smooth decreasing function supported on [0, 1] with Jj^ (j^ix) dx = 
1. Put 



and 



so that 

(4.2) ±^{x) <l{x <0} <i^,. 
Note also that 

(4.3) increases with e for all x 
and 

(4.4) 4'e{x) decreases with e for all x. 

Let b^ and 6^ be the two barriers corresponding to i}} Jyx) and V'g- The existence and 
uniqueness of these two barriers follows by Theorem 12. II From (14.21) we have that 

Ut)<Ut) 

for all t and from (g^) that 

b^{t) is increasing in e for each t 

and 

&g(t) is decreasing in e for each t. 



Put 



6,(t) = lim5,(t) 
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and 



Then, 

(4.5) Kit) < hit), 

and both of these functions give a stopping time with the correct distribution for 
the case where ip is the indicator of (—00, 0]. Because of (|4.5p . it must be the case 
that b^it) = fe^(t) for Lebesgue almost all t. □ 

5. Pricing Claims 
Suppose that the asset price iXt)t>Q is a geometric Brownian motion given by 

(5.1) 



= lidt + adWt, 



where iWt)t>o is a standard Brownian motion. We model default using a diffusion 
iYt)t>o where 

(5.2) dYt ^ dBt, 

with (-Bt)t>o another standard Brownian motion. We assume that the Brownian 
motions W and B are correlated with correlation — 1 < p < 1; that is, the cross- 
variation of the two processes satisfies 

[B,W]t = pt, t>0. 

We can assume without loss of generality that for two independent Brownian mo- 
tions B',B" we have 

f Wt = B'„ 



[ Bt = pB', + v/1 - P'B','. 
In the following we will look at pricing contingent claims with a fixed maturity 
T > and payoff of the form 

F(XT)l{r > T} 

for the random time 

r inf |t > : ^(F^ - 6(s)) ds > C/| , 

where U is an independent exponentially distributed random variable with mean 
one. 

Note that 

ijiY, - bis)) ds 



E''[FiXT)l{T > T}] =E^ 



F(XT)exp ^-^^ 



More generally, we will be interested in expressions of the form 



FiXT)l{T>T} iX,)o<s<t,T>t 



FiXr) exp -A / 7/i(n - bis)) ds (X,)o<.<t, r > t 



which we interpret as the price of the payoff at time < t < T given that default 
has not yet occurred. 
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Consider the Markov process Z — (X, Y, V) where X, Y are as above, and F is a 
process that, when started at is at t; + f after t units of time, that is, Vt = Vo+t. 
The generator of Z is 

A = {l/2)a^x^Dl + nxD^ + (l/2)£»2 + paxD^Dy + Dy. 

We want to compute 



F{XT)e-^o \HYs-b(v,))ds 



The Feynman-Kac formula says that the solution to the PDE 
(5.3) 



Dtu{x, y, V, t) = Au{x, y, v, t) - Xip{y - b{v))u{x, y, v, t), 
u{x,y,v,Q) = F{x), 



satisfies 



E(^'3') F(XT)exp^-^ Xip{Y, - b{s))ds 



u{x,y,0,T). 



Thus, if wc assume the Brownian motion Y has an random starting point Yq with 
density / that is independent of (It — lb)t>o, then 



F{Xt) exp 



X^{Ys - b{s))ds 



u{x,y,0,T)f{y)dy. 



Using this and the Markov property, one can find the function K{x, y, t) satisfying 



K{Xt,Yt,t)=E' 



F{XT)exp l-xj^ ^(n - b{s))ds\ I (X,)o<.<t, (n)o<a<t, r > t 



The price at time t, given that we know the history of the price process Xt and 
that default has not happened up to time t, is 



E 



F{Xt)1{t>T} {Xs)o<s<t, t > t 



E 



K{Xt,Yt,t) iXs)o<s<t, T > t 



E[K{Xt,Yt,t)l{T > t} I {Xs)o<s<t] 
E[l{r > t} I (X,) 

0<s<tj 



It follows from the SDE for X that 

1 " 



B' = W,= 



logXt-logXo+ ( y-M)* 



so if we observe the asset price X, then we can reconstruct the standard Brownian 
motion B' . On the other hand, 



Xt = Xoex^{aB[- — - /x U 
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Now, 

E[K{Xt,Yt,t)l{T>t} I {X,)o<s<t] 

= E 



K i^Xo exp l^aB't - y - /^j ^ j , + pB't + ^1 - P^B'^', t 
1 1^ (Fo + + - P'S;' - 6(s)) < C/ 



We therefore want to be able to compute for a function c 
expected value 



the conditional 



E 



K (xoexp(ac{t)- (^-pUj ,Yo + pc{t) + ^1 - p^B'^' ,t 



(Yo + pc{s) + - p'B" - Ks)) < f^J 
X (^Xo exp (^ac(t) -(^-p^l^^Y^ + pc{t) + ^1 - P^B^ 
X exp ^ ^ (i"o + pc(s) + - P'-B" - c^s^ 



with {B^)t>[) a standard Brownian motion independent of Xq. We can do this 
using Feynman-Kac. 

The denominator in the formula for the price at time t is a special case of the 
numerator we have just calculated with K = 1, and it can be dealt with in the 
same way. 

We have thus observed that computing the price of a contingent claim reduces 
to solving certain PDEs with coefficients depending on the path of the asset price. 



6. Numerical Results 

In this section we present the results of some numerical experiments. We solved 
the PDE/ODE system p.l2p using the pseudo-spectral Implicit- Explicit Fourth 
Order Runge-Kutta scheme ARK4(3)6L[2]SA-ERK, taking 8192 nodes and period 
16, developed in |KC03| . For the function ^jj we used the Fejer kernel of order 
512 applied to the indicator of the set {a; G R : x < 0}; in other words -0 is the 
Cesaro sum of the truncated Fourier series of order 512 of the indicator of the set 
{x S M : X < 0}. The time horizon was taken to be T = 8, the initial distribution 
of the credit index process Y was taken to be normal (Yq ~ (0, a^) with standard 
deviations a = 0.0625,0.125,0.25,0.5), and the time to default was taken to have 
an exponential distribution {G{t) = e""^* with rates v = 0.0625,0.125,0.25,0.5). 
We show the resulting barrier b in Figure [1] We also show the relative error be- 
tween the survival function G{t) and the numerically computed value of J^u{x, t) dx 
(recall (|2.3p ). and the relative error between the hazard rate —g{t)/G{t) and the 
numerically computed value of /jj "0(2; — b{t))u{x, t) dx/ u{x, t) dx (recall (|2.4p ). 
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Figure 1. This figure displays the results of the numerical ex- 
periments described in Section El Successive rows represent the 
values a = 0.0625,0.125,0.25,0.5 of the standard deviation for 
the initial distribution of the credit index process Y . The first 
panel in each row gives the barriers for the rate parameters v — 
0.0625, 0.125, 0.25, 0.5 of the exponential default time distribution. 
The second (resp. third) panel in each row gives the relative errors 
between the actual survival function values G{t) (resp. the actual 
hazard function values —g{t)/G{t)) and the numerically computed 
ones - see the text for details. 



7. Calibrating the default distribution using CDS rates 

For the sake of completeness, we review briefly the scheme proposed in |DPllj 
for determining the distribution of the time to default. 
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A credit default swap (CDS) is a contract between two parties. The buyer of 
the swap makes a number of predetermined payments until the moment of default. 
The seller is liable to pay the unrecovered value of the underlying bond in the event 
of a default before maturity. Normalizing the notional value of the bond to 1, the 
seller's contingent payment is 1 — i?, where R £ (0, 1) is the recovery rate, which 
we take to be constant. The premium payments are made at a set of times {ti}. 
The maturities are a subset of the premium payment times; that is, they are of the 
form To = 0, Tj = tf^f^j-^ , j = 1, . . . , n. For j — 1, . . . ,n there is an upfront premium 
tt'j and a running premium rate nj (having accrual factors 6i). Denote the price 
at time zero of a zero coupon risk- free bond with maturity tj by po{tj). It follows 
from standard non-arbitrage arguments that 



(7.1) 



E 



S,po{U)G{t,) - (1-i?) 



E 

=k{]-i)+i 



where G{t) — P{r > t} is the tail of the distribution of the time to default. 

Suppose now that the default distribution has piecewise constant hazard rate; 
that is, that 

G{t) = exp ^- h{s) ds^ , t> 0, 

where h(s) — hi for s G [Ti,Ti^i). Given the market data (ttijItI), (7r2,7r2), . . . we 
can find, using equation (j7.ip . the constants ho, hi, ... . 

We use the following procedure to find the barrier b. Set v — ho and G{t) = e~^^ . 
Given the initial density /, which we can choose to be any strictly positive function 
/ that is twice continuously differentiable with bounded /, /' and /", we want to 
find a barrier such that for < t < T = Ti we have 

/■* 



/(a;)exp|^— A j ^{x + — b{s)) ds j dx 
This can be achieved by solving the ODE/PDE system (|2.12|) . Next, set i/i = hi. 



T = n 



Ti, fi{x)=E[fix)e^p (-A/o^^V'( 
with 6i(0) = b{Ti) such that on < t < T = 



X -f 

T2 



Bs - b{s))ds 
— Ti we have 



and find a barrier 



fi{x) exp —A / ^^{x + Bg — bi{s)) ds I dx 



This procedure can be repeated until we find a function b on [0, T„] that is con- 
tinuously differentiable everywhere, except perhaps the finite number of points 
Ti, . . . ,Tn. 

8. Duhamel's formula 

For the sake of reference, we provide a statement of Duhamel's formula. Given 
functions v : R x M+ — > R and b : M+ — > R, the solution of 

dx 

I I C/f - 

[ u{x,0)^f{x) xeR, 

is given by 

u{x,t) = [G^{f5t=o)]ix,t) + [G^{~^Pbv)]ix,t) 



dt — ] u — —i^bv, X e R, t > 0, 



KILLED BROWNIAN MOTION 



27 



(8.2) = / G{x-y,t)f{y)dy- / / G{x - y,t- s)tpb(s){y)v{y,s)dyds, 

Jr Jo Jr 

where 

1 »2 

(8.3) G{x,t) := -=e~^, xGR,t>0. 
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